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Abstract — Machine-to-Machine (M2M) communication is now 
playing a market-changing role in a wide range of business world. 
However, in event-driven M2M communications, a large number 
of devices activate within a short period of time, which in turn 
causes high radio congestions and severe access delay. To address 
this issue, we propose a Fast Adaptive S-ALOHA (FASA) scheme 
for M2M communication systems with bursty traffic. The statistics 
of consecutive idle and collision slots, rather than the observation 
in a single slot, are used in FASA to accelerate the tracking 
process of network status. Furthermore, the fast convergence 
property of FASA is guaranteed by using drift analysis. Simulation 
results demonstrate that the proposed FASA scheme achieves 
near-optimal performance in reducing access delay, which outper- 
forms that of traditional additive schemes such as PB-ALOHA. 
Moreover, compared to multiplicative schemes, FASA shows its 
robustness even under heavy traffic load in addition to better delay 
performance. 

I. Introduction 

Machine-to-Machine (M2M) communication or Machine- 
Type Communication (MTC) is expected to be one of the major 
driving forces of cellular networks, as its demand is increasing 
greatly in recent years [1,2]. Behind the proliferation of M2M 
communication, the congestion problems in M2M communica- 
tion become a big concern. The reason is that the device density 
of M2M communication is much higher than that in traditional 
Human-to-Human (H2H) communication [1]. What's worse, in 
event-driven M2M applications, many devices may be triggered 
almost simultaneously and attempt to access the base station 
(BS) through the Random Access Channel (RACH) [3]. Such 
high burstiness can result in congestion and increase response 
time, which motivates our research. 

In literature, several strategies have been proposed for avoid- 
ing radio congestions in M2M communication [4-8]. Among 
these strategies, Slotted-ALOHA (S-ALOHA) type policies, 
e.g., access class barring (ACB) based schemes in [6,7], are 
applied in 3GPP for random access control of M2M devices. 
In these schemes, it is left to users to decide the operation 
parameters such as transmission probability to stabilize and 
optimize the system. Two typical classes of schemes, additive 
and multiplicative schemes, have been proposed for stabiliz- 
ing the S-ALOHA system. Historical outcomes are applied 
in these schemes to estimate the network status and optimize 
the access probability. However, as discussed in more detail 
later, traditional additive schemes, such as Pseudo Bayesian 
ALOHA (PB-ALOHA) [9], estimate the number of backlogged 
devices based on the observation of the previous slot and cannot 
adjust the transmission probability in time under highly busty 
traffic, which results in large access delay. The delay can be 



shortened in multiplicative schemes [10], e.g., Q-Algorithm 
[11] and its enhanced version Q+ -Algorithm [12]. The reason 
is that the estimate has exponential increment in consecutive 
collision slots and exponential decrement in idle slots, which 
means multiplicative schemes can track the network status in a 
short period. However, the throughput suffers in these schemes 
due to the fluctuations in the estimation [10]. More recently. 
Adaptive Traffic Load Slotted Multiple Access Collision Avoid- 
ance (ATL S-MACA) mechanism in [8] uses packet sensing 
and adaptive method to improve the access performance under 
heavy traffic load. But the scheme is designed for M2M com- 
munications with Poisson traffic and is not suitable for event- 
driven M2M applications. 

In this paper, we propose a Fast Adaptive S-ALOHA (FASA) 
scheme for access control of event-driven M2M communica- 
tions. In order to deal with congestions resulted from the bursti- 
ness, we collect access results in the past slots, in particular, 
consecutive idles or collisions, and apply them to track the 
network status. Furthermore, using drift analysis, we carefully 
design the parameters in the scheme such that the transmis- 
sion probability can converge quickly to the optimal value. 
With numerical simulations, we show that the proposed FASA 
scheme can achieve the near-optimal performance in reducing 
the access delay, as well as robust performance under all traffic 
loads less than e^^, which is the maximum throughput of S- 
ALOHA system. 

The remainder of the paper is organized as follows. In Sec- 
tion II, we present the system model, including the bursty traffic 
model for event-driven M2M communications. In Section III, 
after analyzing the limit of traditional fixed step-size adaptive 
policies, we propose the FASA scheme and design it based on 
drift analysis. In Section IV, simulation results are presented to 
evaluate the performance of the proposed scheme, compared 
with the theoretical optimal scheme, PB-ALOHA, and Q+- 
Algorithm. In Section V, we discuss future work and conclude 
the paper 

II. System Model 

In this paper, we consider S-ALOHA based random access 
control for bursty M2M communication traffic. The system 
consists of a BS and a large number of M2M devices, where 
devices with data attempt to access the BS through a single 
RACH. The time is divided into time slots, each of which 
is long enough to transmit a request packet. Deferred first 
transmission (DFT) mode [9] is assumed, in which a device 
with a new request packet immediately goes to backlogged 
state. In slot t, all backlogged devices transmit packets with 



probability pt, which is broadcasted by the BS at the beginning 
of the slot. Moreover, an ideal collision channel is assumed, 
where the transmitted packet will be successfully received by 
the BS when no other packets are being transmitted in the same 
slot. Let Zf denote the access result in slot t, and Zt = 0, 1, 
or c depending on whether zero, one, or more than one request 
packets are transmitted on RACH. At the end of slot t, the BS 
decides the transmission probability for next slot based on the 
sequence {Zq, Zi, . . . , Zt}, i.e.. 



pt+i = nt(Zo, Zi,..., Zt). 



(1) 



The objective of the BS is to maximize the throughput and 
minimize the access delay. It is well known that, when A^t > 1 
in slot t, where Nt is the number of backlogged devices in slot 
t, using transmission probability pt = maximizes the 

throughput of the S-ALOHA system. However, the BS does not 
know Nt and has to obtain its estimate Nt based on the access 
results in the past. 

For the traffic model, we focus on event-driven M2M ap- 
plications, i.e., when an event is detected, a random and large 
number of M2M devices become active almost simultaneously 
and attempt to access the BS on the RACH. We call this an 
active stage. To capture the burstiness of M2M traffic, instead 
of assuming Poisson arrival process, we assume that when an 
event is detected, N devices are triggered in a short duration T. 
As suggested in [6], the active time follows the beta distribu- 
tion, of which the probabiUty density function is 



/(^) = 



13-1 



e [0,r], 



(2) 



of devices attempting to access the channel. To stabilize and 
optimize the S-ALOHA system, Nt needs to track the actual 
number of backlogged devices Nt, especially when Nt is large. 
When Nt ~ n and Nt ~ n, the estimation drift can be 
calculated as [13] 



where B{a,l3) is the beta function, and [6] suggests that the 
values Q = 3 and /? = 4 give the best fit. The number of active 
devices N and the active duration T are both random variables 
and no prior knowledge about them is assumed at the BS. For 
the sake of simplicity, a backlogged device will not generate 
any new requests since the new coming data can be transmitted 
as long as the device accesses the BS successfully. On the other 
hand, all the backlogged devices will keep retransmitting until 
their request packets are successfully received by the BS. 

III. Fast Adaptive S-ALOHA 

The estimation of the number of backlogged devices plays 
an important part in stabilizing and optimizing the S-ALOHA 
system. In this section, using drift analysis, we first examine the 
limit of traditional fixed step-size estimation schemes. Then, we 
propose and analyze a fast adaptive scheme, referred to as Fast 
Adaptive S-ALOHA. 

A. Drift Analysis of Fixed Step-size Estimation Schemes 

Many additive schemes with fixed step-size have been pro- 
posed to estimate the number of backlogged devices. A unified 
framework of these schemes is proposed and studied by Kelly 
in [13], where the estimate Nt is updated by the recursion 

Nt+i = m&x{l, Nt+aoI[Zt = 0]+aiI[Zt = l]+aJ[Zt = c]}, 

(3) 

where qq, ai, and Gc are constants and I[A] is the indicator 
function of event A. 

With the estimation, the BS sets the transmission probability 
Pt = l/Nt for all backlogged devices, and thus the offered 
load p — NtPt = Nt/Nt, representing the average number 
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as n ^ cxD, with njh — p fixed. 

By properly choosing the parameters (i = 0, l,c) such 
that A(p) < if p < 1 and A(/9) > if p > 1, the 
estimate Nt will drift towards the true value and thus the S- 
ALOHA system can be stabilized. However, these fixed step- 
size schemes are not suitable for systems with bursty traffic. 
When the estimate Nt deviates far away from the true value Nt, 
we have limp_).o = ao and limp^oo ^(p) = Oc- These 

limits indicate that the drift tends to be a constant even when 
the deviation is large, which could result in a large tracking 
time. Thus, it is necessary to design fast estimation schemes 
for event-driven M2M communication. 

B. Design of FASA 

As analyzed in the previous subsection, fixed step-size es- 
timation schemes such as PB-ALOHA may not be sufficient 
to adapt in a timely manner for systems with bursty traffic 
because it uses a constant step-size even when the estimate is 
far away from the true value. We note that in addition to the 
access result in the previous slot, the access results in several 
consecutive slots will be helpful for improving the estimation 
as they may reveal additional information about the true value. 
Intuitively, collisions in several consecutive slots are likely 
caused by a significant underestimation, i.e., Nt <C Nt, and the 
BS should aggressively increase its estimate. Similarly, several 
consecutive idle slots may indicate that the estimate Nt ^ Nt, 
and it should be reduced aggressively. 

Motivated by this intuition, we propose a FASA scheme that 
updates Nt as follows 

max{l, Sft-l- hoii^XKo^tY}, if Zt = 

Nt^, = i if ^* = 1 



Nt 



1 



+ h,{iy)iK,,tr, 



if Zt 



(5) 

where K^ t and K^.t are the numbers of consecutive idle and 
collision slots up to slot t, respectively; > is the parameter 
that controls the adjusting speed; ho{iy) and hc{iy) are functions 
of 1/ that guarantee the tracking criterions [13]. We are mostly 
interested in those cases where the number of backlogged 
devices is large. Hence we will approximate max{l,a-} as x 
for notational simplicity in the analysis later. 

Next, we design ho{i^) and hc{i^) by drift analysis. Consider 
tracking a fixed number of backlogged devices, i.e., Nt = n is 
constant for all t £ Z_|_ {0, 1,2,.. .}. Assume that in slot t, 
the estimate Nt = n, and thus the offered load p = n/n. When 
n is large, the drift of estimation for FASA can be calculated as 

Afasa(p) = E[Nt+i - Nt] 
= qo[p)E[/^^{p)]+q^{p)E[^M+1c{p)E[^,{p)l (6) 



where qoip) ^ e p, qi{p) = pe p, and qdp) = 1 - qo{p) - 
qi [p) are the probabiHties of an idle, success, and collision slot, 
respectively; Ai{p) (i = 0, l,c) is the change in Nt resulting 
from the corresponding update. 

Obviously, Ai{p) = since the estimated number remains 
unchanged when a packet is successfully transmitted in a slot. 
However, Ao{p) and Ac(p) depend on the distribution of KQ,t 
and Kc,t, respectively. For example, 

t+i 

E[Ao{p)] = E[Aoip)\Ko,t = ko]PriKo,t = h). 

Unfortunately, it is difficult to obtain the distribution of 
Ko,t and Kc,t directly, and thus we resort to approximation 
in order to make the problem tractable. It is easy to see that 
the update of Nt is mostly affected by the access results in the 
past s slots {Zt^s, Zt-s+i, Zt-i}, where s is the number 
of consecutive idles or collisions immediately proceeding t, 
or s = 1 if Zt_i = 1. The update step is reset whenever 
Zt-i 7^ Z t oi Zt = 1. For tractability, we approximate the 
estimates Nf {t' = 0, 1, . . . , t — 1) as Nf Thus, the access 
outcomes {Zq, Zi, Zt} are independent and identically dis- 
tributed (i.i.d.) random variables, of which qa{p), qi{p), and 
qdp) are the probabilities of the idle, success, and collision slot, 
respectively. As we will see later, this provides a rough bound 
for the drift and guarantees the convergence of the algorithm. 

First, to calculate the expected drift in an idle slot, suppose 
that no packet is transmitted in slot t. Then the estimated 
number will be reduced by 1 + hQ{i'){Ko,tY . Ko t = kg 
(ko < t) holds when slots t — fco + 1, t — fco + 2, . . . , t — 1 
are all idle while slot t — ko is not. Under the approximation 
that the access results are i.i.d. in each slot, and the probability 
that an idle slot occurs is qo{p), we have 



Pr{K 



ko) 



q'o"-\p)[l - qoip)], 
cio"'\p), 



if l<ko<t 
if ko = t + l. 



As t tends to infinity, qo{p) tends to 0, and Ko.t can be 
approximated by a geometrically distributed random variable 
with success probability 1 — qo{p)- Hence, 



E[Aoip)] « E - hoi>y)k^o]'lo"'\p)[l qoip)] 



ko = l 



= -[l + ho{l^)M{l^,qo{p))l 



(7) 



where M{v,qo{p)) is the approximate expectation of [Ko.t) 
when the offered load is p, which is given by 



M{v,qo{p)) = J2 k^o<lo°'\p)[l - qoip)] 



(8) 



fco = l 



Secondly, we can calculate the drift of the estimation in a 
collision slot in a similar fashion as follows: 



E[A,ip)] « (e - 2)-i + h,{i^)M{i^,q,ip)). 



(9) 



From (8) and (9), the drift of FASA can be approximated 
using 

ApASAip) w -go(/o)[l + hoiu)M{u,qoip))] 

+qcip)[ie - 2)-l + h,{iy)M{^, q,{p))]. (10) 



In order to keep the estimated number Nt close to the true 
value, it is required that ApASAip) = for p = 1. Letting 
9o = 90(1) = and q* = qdl) = 1 - 2e~^, we expect that 

ApASAil) - ~hoiv)q^Miv,q*o) + hdv)q*aM{v,ql) = 0. 

(11) 

Therefore, in order to satisfy the condition in (1 1), we can select 
the following hoiv) and hdv): 



ho{u)^^l[q*oMiv,q*o)]-\ 

hdv)^'n[q;M{v,q:)]-\ 



(12) 
(13) 



where rj > is a constant and is another parameter for 
controlling the tracking speed. 

We now interpret the effect of i.i.d. approximation. We focus 
on the cases of p <C 1 and p 1, since the evolution of Nt 
in these cases diverges the most from the i.i.d. approximation. 
Up = n/n <^ I in slot t, it is more likely that there were 
consecutive idle slots and the estimate Nt was decreasing in 
the past slots, suggesting Nf > Nt for slot t' < t with 
high probability. In this case, the idle probabilities in the past 
slots would have been larger than that in slot t, while the 
collision probabilities would have been smaller than that in 
slot t. Hence, using Nt as an approximation of Nt' it' < t) 
overestimates E[Ao{p)] and E[Acip)], and the drift analysis 
above gives an upper bound of the drift for p ^ 1. Similarly, 
when p ^ 1, the drift analysis above gives a lower bound 
of the drift. Subsequently, the above analysis using the i.i.d. 
approximation roughly bounds the evolution of Nt in both 
directions and guarantees the convergence of Nt to n following 
Proposition 1 (Section 111-C). 

C. Drift Analysis of FASA 

The chosen hoii^) and hdi^) guarantee that Apasa{^) = 
and thus provide a necessary condition for FASA to track the 
number of backlogged devices. Furthermore, Proposition 1, 
which can be proved by examining the derivative of A fas A (p) , 
shows a desirable property of FASA. With this property, the BS 
inclines to decrease Nt when Nt > Nt and increase Nt when 
Nt < Nt- Thus, Nt is able to track the number of backlogged 
devices n. 

Proposition 1 : Given that hoii') and hdi^) are defined in 
(12) and (13), respectively, the approximate drift of FASA 
ApASAip) is a strictly increasing function of p. In addition, 
Afasa{p) < when < p < 1 and ApASAip) > when 
p>l. 

Proof: See Appendix. 

In order to understand better the behavior of the scheme, 
we now present the drift of estimation for = 1,2 and 3. 
When 1/ G Z+, Mil/, qdp)) [i = 0,c) is the i^th-moment 
of a geometrically distributed random variable with success 
probability 1 — qi{p), and its closed-form expression can be 
obtained. Consequently, with the expressions of Mil', qdp)), 
we can obtain functions hoii'), hdv), and the drift ApASAip) 
for = 1, 2 and 3. For example, when v = 2, the drift of FASA 
can be expressed as 



A 



(2) 

FASA 



(p) = A(")+7y 



(e-l)2(eP + l) 



ie + l){eP - 1)2 
2(e''-p-l)(2e''-p-l 



(e-2)(e-l)(p+l)2 



(14) 



where A*^*^^ ~ —e p + ^ is the drift of the fixed step-size 
estimation (3) with oq = —1, ai = and Oc = (e — 2)~^. 

As shown in Fig. 1, when the estimated number deviates far 
away from the actual number of backlogged devices, FASA 
adjusts its step-size accordingly, while PB-ALOHA still uses 
the same step-size. Therefore, using FASA results in much 
shorter adjusting time than PB-ALOHA, and thus improves 
the performance of M2M communication systems with bursty 
traffic. Note that the drifts of multiplicative schemes such as 
-Algorithm are not illustrated here since they depend on not 
only the offered load p but also the estimate Nf 
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Fig. 1. Drift of estimation (r? = 1) 

IV. Simulation Results 

In this section we evaluate the performance of the proposed 
scheme through simulation. We compare the performance of 
our FASA scheme, the ideal policy with perfect knowledge of 
backlog, PB-ALOHA [9], and Q+ -Algorithm [12]. 

A. Settings 

First, we choose the parameter values for the adaptive 
schemes. Since simulations show that the performances of 
FASA with different rj and v are close, we just present the 
results for 77 = 1 and v = 2 \n this paper due to the space 
limitation. With perfect knowledge of Nt, the ideal policy sets 
transmission probability at pt = 1 /Nt for Nt > 0. Thus, the 
ideal policy achieves the minimum access delay of S-ALOHA 
and serves as a benchmark in the comparison. For PB-ALOHA, 
we use the estimated arrival rate At = e^^, as suggested in [14]. 
Q"*" -Algorithm belongs to the class of multiplicative schemes 
which is first proposed by Hajek and van Loon [10]. In Q+- 
Algorithm, Nt is updated as follows: 

Nt+i = max{l, [I{Zt = 0)/Co+I{Zt = l)+CcI{Zt = c)]Nt}, 

where Co = '2°-^^ ~ 1-1892 and Q = 20-35 ^ ;l 2746 are 
suggested in [12] for optimal performance. 

Second, we describe the simulation scenarios. In order to 
gain more insights into the performance of access schemes with 
bursty traffic, we present simulation results for a single active 
stage and assume that the intervals between two consecutive 
active stages are sufficiently large so that all the devices can 
successfully access the BS before the next trigger. During the 
active stage, devices are triggered according to the beta 
distribution with parameters a = 3, /3 = 4 [6], and time span 



T = 50 slots, which is the number of slots in one second when 
the PRACH period [3] is 20 ms. It is assumed in 3GPP that 
30,000 devices can be triggered in 10 s [6], thus we choose N 
to be 100 to 3,000 per second. 

B. Results 

In some event-driven M2M applications, response can be 
taken with partial messages from the detecting devices and not 
all devices need to report an event. Thus, both the distribution 
of access delay and average delay are evaluated to study the 
performance of the proposed scheme. 

Fig. 2 shows the cumulative distribution function of access 
delay for different schemes. From this figure we can see that 
the performance of the proposed FASA scheme is close to the 
benchmark with perfect knowledge. For PB-ALOHA, it takes 
a long time to track the number of backlogged devices and few 
devices can access successfully during this period. For instance, 
the 10% delay, which is the access delay achieved by 10% 
of the active devices, is much larger than other schemes. For 
example, when = 500, the 10% delay of FASA is about 280 
slots while it is 520 slots for PB-ALOHA. With multiplicative 
increment, the Q+ -Algorithm can track the number of backlogs 
in a short time because of the exponential increment due to the 
consecutive collision slots. However, it takes longer for all the 
devices to access the channel under Q+ -Algorithm than under 
FASA due to the large estimation fluctuation in Q^-Algorithm. 
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Fig. 2. Cumulative distribution function of access delay 

Fig. 3 compares the average delay of the access schemes. 
The top figure indicates that the average delay increases almost 
linearly as the number of active devices increases. To quantify 
the divergence from the theoretical optimum performance, we 
define the normalized divergence as e{D) = ^jj? ^ 100%, 
where D is the average delay of a particular scheme and D* 
is the theoretical optimum delay with perfect information. The 
bottom part of Fig. 3 shows the divergence of PB-ALOHA, 
Q+-Algorithm, and FASA. The delay of PB-ALOHA scheme is 
larger than the optimum value by about 22%. The divergences 
of Q+ -Algorithm and FASA are both much less than that of 
PB-ALOHA. As the number of the active devices increases, 
the effect of estimation fluctuation becomes small and the 
performance gets close to the optimum value. For large N, the 
divergence is about 5% for Q+ -Algorithm and 2% for FASA, 
which indicates that FASA scheme performs slightly better than 
Q+ -Algorithm. 



As shown in Figs. 2 and 3, for a single active stage, the 
performance of FAS A and Q+ -Algorithm is close, with the 
proposed FASA scheme performing slightly better However, 
when considering the long term performance, as discussed in 
[10], for given Co and ^i, the estimated value Nt continues to 
fluctuate when it gets close to Nt and the stable throughput is 
less than some maximum value 8max < e^- We verify this 
property through simulations for Poisson arrival process with 
mean A. The results show that, for PB-ALOHA and FASA, 
the number of backlogged devices is finite though the average 
value becomes large as A becomes close to e~^. For Q+- 
Algorithm with the given parameter values, however, when A 
is larger than about 0.36, the number of backlogged devices 
grows unbounded, indicating that the algorithm is unstable 
when A > 0.36. 

Average delay 



and 

9cip) = \qcip)E[A,{p)]\ = q,{p)[- 
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Fig. 3. Average delays for different schemes 



V. Conclusion and Future Work 

This paper proposes a fast adaptive S-ALOHA scheme, 
called FASA, for event-driven M2M communications. An ap- 
proximate drift analysis and simulation results show that using 
FASA, a BS can track the number of backlogged devices more 
quickly. That is a main advantage compared to fixed step-size 
additive schemes, e.g., PB-ALOHA. Compared to multiplica- 
tive schemes, simulation results show that the proposed FASA 
scheme has better stability performance under heavy load in 
addition to slightly better delay performance. Currently, the 
analysis presented in this paper is based on the approximation 
with fixed number of backlogged devices. A rigorous analysis 
about the properties of the proposed scheme with general arrival 
traffic, e.g., mixed bursty and Poisson traffic, is left as our future 
work. 
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Appendix: Proof of Proposition 1 
The proposition can be proved by calculating the derivative 

of ApASAip)- 

For a given value of i/, let 

9oip) = \qo{p)E[Ao{p)] \ ^ qoip)[l + hoiiy)M{iy,qoip))], 



hc{u)M{u, qc{p))]. 



Thengo(l) = 3c(l) = e ^ +77 and Afasa(I) = -3o(l) + 
.gc(l) = 0. Next, we claim that, for given 1/ > 0, M{v^ q) = 
SfeLi fc'''J''^~^(l — 5) is an increasing function of g (0 < q<l), 
because 



dM 



fe=i 

k* 

Y,k'^q^-^[k{l-q)-l] 



CXD 

00 

> {k*YY^q''-^[k{\-q)-l]=0, (15) 
fc=l 

where k* ~ \_T^\ is the largest integer not greater than j-^^, 
and thus k" < {k*Y if 1 < fc < fc* and fc"^ > {k*Y 
if k > fc*. In addition, the idle probability qa{p) = e^P is 
nonnegative and strictly decreasing in p. Hence, M{i', qo{p)) 
is strictly decreasing in p and ,90 (p) is a strictly decreasing 
function of p. On the other hand, since qdp) = l — e~P~pe~P 
is nonnegative and strictly increasing in p, we can show sim- 
ilarly that gc{p) is a strictly increasing function of p. Thus, 
Afasa{p) = —go{p) + gc{p) is a strictly increasing function 
of p. Consequently, Afasa{p) < Apasa{^) = when 
< /9 < 1 and Afasa{p) > ApASAi'^-) = when p > 1. 
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